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Abstract. For the two dimensional kinetic Ising model at finite temperature, the 

local mean magnetisation Mt = t~ 1 I a(t')dt' , simply related to the fraction of time 

J 

spent by a given spin in the positive direction, has a limiting distribution, singular at 
±mo(T), the Onsager spontaneous magnetization. The exponent of this singularity 
defines the persistence exponent 9. We also study first passage exponents associated 
to persistent large deviations of Mt, and their temperature dependence. 



submitted for publication to J. Phys. A: Math. Gen. 

PACS numbers: 02.50.Ey, 05.40. +j, 05.50+q, 75.10.Hk 

In this work we present a new approach to the study of persistence for systems 
undergoing phase ordering [jl], ^ at finite temperature, which we shall illustrate on the 
case of the two dimensional Ising model. In this approach persistence appears as a 
stationary property of the coarsening system, and the role of the Onsager spontaneous 
magnetization at equilibrium mo(T) is made apparent, thus revealing new fundamental 
features of phase ordering. It departs from previous approaches to finite temperature 
persistence [|[ |[ |[ [|, where these features did not appear. 

Consider a system of Ising spins o"j(f) = ±1 located at sites i = 1, . . . ,N, started 
from a random initial condition, and evolving under the heat bath dynamics at fixed 
temperature below the critical temperature. At each time step a spin is picked at 
random, and updated with the probability 



/"'(-,!/ -!-(!/! = +1) - | 1+tanh-^a (/) j (1) 

where the sum runs over the neighbours of site i. Under this dynamics spins thcrmalizc 
in their local environment. Therefore the system coarsens, i.e. domains of opposite 
signs grow and, in the scaling regime, the system is statistically self similar, with only 
one single characteristic length scale, which is the size of a typical domain jl], ||] . 

The question of persistence is to determine the fraction of space R(t) which 
remained in the same phase up to time t [j?], ^ (or from time t\ to time £2)- For the 
two dimensional Ising model at zero temperature, two phases coexist, corresponding 
to all spins equal to +1 or all spins equal to —1. Hence R(t) is equivalently defined as 
the fraction of spins which did not flip up to time t M, i.e. which were not swept by 
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an interface between domains of all spins +1 or all spins — 1. Numerical measurements 
indicate an algebraic decay R(t) ~ t~\ with 6 « 0.22 || ^ [□]]. 

The definition of persistence at finite temperature below T c is more subtle to 
implement because one has to make clear what is meant by 'phase'. By essence, in 
the coarsening process there is phase separation, each phase wanting to develop at the 
expense of the other. It is therefore intuitive that the system, though perpetually out 
of equilibrium, tries to reach locally one of the two equilibrium phases, corresponding 
to ±?7io(T), where 

2 N ' 



m (T)= [smh T j j (2) 

is the Onsager spontaneous magnetization at equilibrium |12| . Hence in the scaling 
regime the average magnetization inside a domain measured on a scale of time small 
compared to the flipping time of the domain, should be close to the equilibrium 
magnetization at this temperature. Coming back to persistence, the definition of 
R(t) should reflect, in one way or the other, the fact that a given point in space 
remained in a phase of average magnetization equal to ±mo(T) up to time t. This 
intuitive analysis is confirmed by what follows. 

The central point of our approach is to consider the statistics of the local mean 
magnetization -simply related to the fraction of time spent by a spin in the positive 
direction- in the limit of large times. This line of thought was already used in (l^] in the 
study of domain coarsening for the one dimensional Ising model at zero temperature, 
or for the simple diffusion equation evolved from a random initial condition (see also 

0) . The idea is that since persistence probes the past history of the system, a 
natural quantity to consider is f Q a(t')dt' — T ( + — T t ~, where a(t) is the spin at site 
i and T t + (T t _ ) is the length of time spent by the spin pointing upward (downward), 
with t = T t + + TfT . The local mean magnetization is defined as 

M t = \ [ a(t')dt' = 2^- - 1. (3) 
1 Jo t 

For instance at zero temperature, the persistence probability R(t) is equal to V(M t = 

1) since the event {<r(t') = 1,W < t} is identical to the event {M t — 1}. Since M t 
is a local quantity varying from site to site, one is naturally led to investigate the 
distribution of Mt, 

P(t, x) = V(M t > x) {-l<x<l). (4) 

For the one dimensional Ising model at T = is was shown in ]l3| by analytical 
arguments and numerical measurements that, when t — * oo, P(t, x) converges to a 
limit distribution P^a;) with density 

fu{x) = -^Poo(x) (5) 

singular at x = ±1 with singularity exponent equal to 9 — 1. It is for instance easy to 
show that, when t — > oo, the limit of (M t 2 ) is a constant equal to ^4(1), the Laplace 
transform of the two-time correlation with respect to the variable lnt, at argument 
equal to 1 |l3[ ] . This result therefore provides a stationary definition of persistence at 
zero temperature. The same holds for the diffusion equation O, H. 
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Figure 1. Histogram of the density of Mt at time t = 1000, for T = 
0, 0.8T C , 0.9T C , 0.95T C , T c , 1.1T C , from bottom to top (at the central point x = 0). 
(The system size is 3072 2 .) 



We now address the same questions at finite temperature. We first report on 
numerical results. Figure 1 depicts the histogram of the density of M t at time t — 1000 
for values of T ranging from to 1.1T C . Already for such a short time, and for 
every temperature T < T c , the density is maximum around ±mo(T), the equilibrium 
magnetization (||). At larger times and for T > T c , the density of M t becomes peaked 
around zero, i.e. the mean magnetization converges toward the average magnetization 
per spin (a) = 0, reflecting the fact that the system reaches equilibrium. At T c the 
peaking of the density of M t is observed to be very slow. Finally at T < T C) P(t, x) 
converges, when t — ► oo, to a limit distribution Pco(x) with density Jm(x) given as in 




The existence of a limit law at finite temperature T < Tc relies on the same 
arguments as for the zero temperature case. For instance the convergence of (A/ 4 2 ) to a 
constant equal to A(l) still holds since it only relies on the existence of a scaling regime 
p3[ . The striking fact is that now the density concentrates on [— mo(T), mo(T)] with 
an exponential decay with time of V(M t > too(T)) to 0. Moreover the limit density is 
singular at ±mo(T), with a singularity exponent 9 — 1 which defines the persistence 
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Figure 2. Plot of 1 — P(t, x) against the rescaled variable x/mo(T), at T = (full), 
T = 0.8 T c (dots), for t = 20 000, and at T = 0.98 T c , for t = 1000, 5000 and 30 000. 
(The system size is 3072 2 .) 

exponent at temperature T. 

This leads to the question of the temperature dependence of persistence for T < T c . 
The simplest assumption is that, during the coarsening process, the scales of time 
between a short time regime and the scaling regime decouple, yielding the following 
relation between the moments of the limit distributions at T and at zero temperature, 

{M 2k ) T = (m (T)) 2k (M 2k ) (6) 

implying the identity of the limit distributions, if M t is rescaled by mo(T), and as a 
consequence, the temperature independence of 9. This would be in agreement with 
the usual view that zero temperature is an attracting fixed point for the dynamics 
of phase ordering j^, Q|. Equation (^) is hard to check by numerical measurements 
because the convergence of the data is observed to be slow. The difficulty is illustrated 
by figure 2 which depicts a plot of 1 — P{t,x) against the rescaled variable x/mn(T), 
at T = and 0.8 T c , for t = 20 000, and at T = 0.98 T c , for t = 1000, 5000 and 30 000. 
Though one cannot be conclusive on the sole basis of numerical measurements, data 
collapse seems nevertheless plausible. Let us note that the limit distribution 1 — (x) 
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at T = is well approximated by a beta distribution, as was observed for the Id Ising 
model 13 1, or for the diffusion equation jl3[ |l4j]. The singularity exponent of the 
beta distribution is found to be around 0.22. A more precise numerical determination 
of the exponent 9 from the limit distribution of M t needs further work and will be 
presented elsewhere. 

Let us summarize at this point. For T < T c , the local mean magnetization Mt has 
a limit probability density when t — > oo, defined on the interval [— mo(T), mo(T)], 
where mo(T) is the equilibrium magnetization, and singular at both ends. This 
provides a stationary definition of persistence, which is a natural extension of the zero 
temperature case, where the singularity exponent defines the persistence exponent. 
These are the central results of the present work. 




Figure 3. Log-log plot of R(t, x) for the two dimensional Ising model at zero 
temperature, x varying from —1 to 1, by steps of 0.1, from top to bottom. (The 
system size is 1536 2 .) 

Another new aspect of persistence introduced in |l3[ ] is concerned with persistent 
large deviations. The probability of persistent large deviations above the level x 
(— 1 < x < 1), denoted by R(t,x), is defined as the probability that the mean 
magnetization was, for all previous times, greater than x [|l3| , 

R(t,x) = V(M t i > x,W < t). (7) 
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In other words one is interested in the persistence probability of the stochastic process 
a(t, x) = sign(M t 



x) 13 



If one views the stochastic process a(t) as the successive 
steps of a fictitious random walker, then M t is the mean speed of the walker between 
and t, and R(t, x) is the probability that this mean speed remained larger than 
x between and t. This probability is a natural generalization of the persistence 
probability R(t) = R(t, 1), which corresponds for the walker to always stepping to the 
right. 

For the one dimensional Ising model at zero temperature, R(t, x) was observed to 
decay algebraically at large times with an exponent 9{x) continuously varying with x 
Q . For x = 1, 9(1) — 9, the usual persistence exponent. Figure 3 depicts a log-log 
plot of R(t, x) for the two dimensional Ising model at zero temperature, x varying 
from —1 to 1, while figure 4 depicts the corresponding exponents 9(x), extracted from 
figure ||. Let us mention that algebraic decay of R(t, x) was also observed for the 
diffusion equation |l3| , and that this quantity and the corresponding exponents 9(x) 
can be exactly computed for the simple model considered in [191. 




Figure 4. Exponents 9(x) extracted from figure 3. 



We now address the role of temperature for persistent large deviations. As is 
obvious from the first part of this work, if x > mo(T), then R(t,x) decays to zero 
exponentially rapidly. On the other hand, for x < mo(T) one observes algebraic decay 
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of R(t, x), as at zero temperature. Otherwise stated, x = m (T) separates two regimes 
of persistent large deviations, between exponential and algebraic. As a consequence, 
and by analogy with the zero temperature case, one could think of extracting the 
persistence exponent at finite temperature from the decay at large times of R(t, x) 
when x — > mo(T) _ , which leads to the formal definition 



R(t)= lim R(t,x) 



= lim 9(x). 

x-»m (T)- 



(8) 



However this definition is not easy to implement in practice. 




Figure 5. Log-log plot of R(t,0) for temperatures T = 0, 0.8T C , 0.9T C , 
0.95T C , 0.98T C , T c , 1.1T C , 1.5T C , from top to bottom. (The system size is 1536 2 .) 



We did not investigate the temperature dependence of 8(x) in all generality. We 
restricted our study to the case x = which corresponds to 



R{t, 0)=v(J^ (j(t')dt' > 0,W <t\= T{T+ > T t 7, W < t). (9) 
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We find that R(t, 0) 



t~ ^ with 



0(0) 
0(0) 
6(0) 



-0.5 

-0.27 

-0.20 



T 
T 
T 




(10) 



For T c < T, after a crossover, 9(0) takes the high temperature value 0.5 while for 
T < T c it takes the low temperature value rj 0.20. (See figure 5.) The explanation 
of the value of 9(0) for T = oo is simple. Since spins are independent, identifying as 
above cr(t), the spin at site i, to the steps of a fictitious one dimensional symmetric 
random walker, R(t, 0) represents the probability that the walker did not cross the 
origin up to time t, which is, as is well known, decaying as t^ 1 / 2 jl6). For decreasing 
temperatures, spins become more correlated, hence the exponent 9(0) decreases. Note 
that the first passage exponent 9(0) is defined for T > T c , i.e. even in absence of 
coarsening. 

This work raises a number of questions. For instance, what is the temperature 
dependence of the two time correlation in the scaling regime, for T < T c l Is the 
hypothesis (^|) valid? What is the behaviour of the distribution of M t at T c when 
t — > oo? At T c , is 9(0) a new independent critical exponent, or is it related (equal?) to 
the persistence exponent 9 C for the global magnetization jl7[ [|? Let us mention that 
for the 3d Ising model the quantities studied here have similar behaviour. Finally, in 
our view, an important point of the analysis presented here is that it may be applied 
to any coarsening system, since it relies mainly on scaling. 
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